Discrimination of crop varieties spanned over heterogeneous agriculture land is a vital application of polarimetric SAR images for agriculture monitoring and assessment. The covariance matrix of polarimetric SAR images is observed to follow a complex Wishart distribution for major classification tasks. It is true for homogeneous regions, but for heterogeneous regions, the covariance matrix follows a mixture of multiple Wishart distributions. We aim to improve the classification accuracy when the terrain under observation is heterogeneous. For this purpose, Wishart mixture model is employed along with expectation-maximization (EM) algorithm for parameter estimation. Elbow method helps us to devise the number of mixtures. The convergence of the EM algorithm depends on the choice of initial points. So, to improve the robustness of the model, different initialization approaches, such as random, K-means, and global K-means, are embedded in the EM algorithm. Further, the degrees of freedom is one of the crucial parameters of Wishart distribution. Therefore, the impact of different degrees of freedom is analyzed on classification accuracy. The method that is equipped with initialization technique along with optimum degrees of freedom is assessed using three full polarimetric SAR data sets of agriculture lands. The first two are benchmark data sets of Flevoland, Netherlands, region acquired by AIRSAR sensor, and third is our study area of Mysore, India, acquired by RADARSAT-2 sensor.